Interaction of Diatomic Molecules with Photon Angular Momentum 
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The interaction of a diatomic molecule with photons carrying well-defined angular momentum 
and parity is investigated to determine whether photon absorption can induce molecular rotational 
transitions between states having angular momentum AJ > 1. A transformation from laboratory 
coordinates to coordinates with origin at the center-of-mass of the nuclei is used to obtain the 
interaction between the photons and the molecule's center-of-mass, electronic, and rotational degrees 
of freedom. For molecules making transitions between rotational levels, there is a small parameter, 
ka -C 1, where k is the photon wave vector and a is the size of the molecule, which enters into the Ej 
and Mj photon absorption probabilities. For electric photons having arbitrary angular momentum 
jh, the probability of absorbing an E(j + 1) photon divided by the probability of absorbing an 
Ej photon, scales as (ka) 2 /(2j + l) 2 . The probability of absorbing an Mj photon, divided by the 
probability of absorbing and Ej photon scales according to the same factor. 



I. INTRODUCTION 

A molecule is significantly more complicated than an 
atom because, in addition to center of mass motion, a 
molecule has vibrational and rotational degrees of free- 
dom. Work on molecular spectra has a long history. A 
significant milestone was made in 1927 by Born and Op- 
penheimer who assumed that the lighter electrons 
in a molecule adjust adiabatically to the slower motion 
of the heavier nuclei, which remain at all times in their 
instantaneous ground state. Further progress in under- 
standing molecular spectra was made by treating a rotat- 
ing molecule as a symmetric top or asymmetric top [IJ, 
sometimes called a rigid rotor d, 0] , see Ref. @ and ref- 
erences contained therein. In the rigid rotor treatment 
of molecules, vibrational and center-of-mass motion are 
essentially neglected. 

Solving the complete quantum mechanical many body 
problem consisting of electrons and nuclei for a given type 
of molecule will in-principle provide the energy spectrum 
of the molecule. In practice, this is a complicated task 
which usually results in heavy numerical work. For this 
reason, much analytic work has been done on simple di- 
atomic molecules [3j, |J] . For modern reviews on the the- 
ory of the molecular Hamiltonians, see Meyer jg], and for 
an overview about diatomic molecules, see Brown and 
Carrington Q- 

While the spectrum of a molecule can in-principle be 
obtained by solving the corresponding quantum mechani- 
cal many body problem, in experiments we usually probe 
molecules with electromagnetic fields. In order to in- 
terpret experiments, we must calculate selection rules 
and probabilities of transition between molecular quan- 
tum states. This involves investigating the coupling of 
a molecule to the electromagnetic field. The fundamen- 
tal theory of a second quantized electromagnetic radia- 
tion field coupled to charges is discussed in Ref. [|| Q . 
Detailed treatments dealing with radiation coupling to 
molecules are given by Atkins and Wooley 1101 . Wool- 
ley [ll|, [HJ , Craig and Thirunamachandran jl3| | and a 



recent detailed treatment is given by Sindelka [14] . 

In the past two decades, a completely different line 
of investigation has dealt with the angular momentum 
carried by the electromagnetic field. Beth [l5| made 
the first experimental observation that light carries an- 
gular momentum. Beth demonstrated that if a beam 
of right-hand circularly polarized light passed through a 
birefringent plate, and the beam was converted into a 
left-hand polarized beam, the plate experienced a torque 
of 2h for each photon in the beam. This experiment is 
cited as a demonstration that a light beam carries spin 
angular momentum. Years later, in 1992, Allen et al. 
showed that laser light with a Laguerre-Gaussian am- 
plitude distribution has well-defined orbital angular mo- 
mentum [161] . Many theoretical and experimental inves- 
tigations followed, see Ref. [l7|. For paraxial beams, i.e., 
in the paraxial approximation, it has been demonstrated 
that the total angular momentum can be separated into a 
sum of spin angular momentum and orbital angular mo- 
mentum [lq |. However, for a general electromagnetic ra- 
diation field, spin and orbital angular momentum cannot 
be separately defined, and only the total angular momen- 
tum is a well-defined quantity 0, H, ■ The interaction 
of atoms with light carrying angular momentum were ex- 
plored theoretically in a number of works [T9M21I ] . 

Within the rigid rotor model of a molecule, one can eas- 
ily see that an electromagnetic radiation field can induce 
(dipole, quadrupole, etc..) rotational transitions in a di- 
atomic molecule through its multipole moments, see for 
example p. 601 in Section 136 of Davydov[4j. However, 
the rigid rotor treatment neglects the molecule center- 
of-mass degree of freedom. In the absence of a radia- 
tion field, the molecular center-of-mass motion decouples 
from the other degrees of freedom and there is no issue 
to consider further, see Eqs. (|3"0|) and (f3"5|) . However, in 
the presence of a radiation field, the molecular center of 
mass plays a key role in the molecule-field interaction, 
see Eq. |g2J) and Ref. dHH^. References dis- 
cuss the non-trivial role of the center of molecular mass 
when a classical plane polarized electromagnetic wave 
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(corresponding to photons having j = h) interacts with a 
molecule. In order to investigate in detail molecular ab- 
sorption of photons having arbitrary angular momentum 
(J > nh), where n > 1, one has to revisit the inter- 
action Hamiltonian between the molecule and radiation 
field. As already remarked, this is complicated because a 
molecule has several degrees of freedom: center-of-mass 
motion, electronic motion, vibrational motion, and rota- 
tional motion. 

In this manuscript, I consider a simple diatomic 
molecule interacting with photons carrying arbitrary an- 
gular momentum, jh. I investigate the mechanism that 
leads to transitions between molecular rotational levels, 
with emphasis on the role of the molecular center-of- 
mass. Of particular interest are probabilities of tran- 
sitions between rotational levels whose angular momen- 
tum differs by more than one unit, A J > 1. I write 
down a Hamiltonian for a diatomic molecule, for the ra- 
diation field, and for their interaction. By making suc- 
cessive transformations from the laboratory coordinates 
into the center-of-mass coordinates, and then to center- 
of-molecule coordinates, I obtain interaction terms that 
show the coupling of a diatomic molecule to the radia- 
tion field. For the molecule, there is a small parameter, 
ka <C 1, where k is the wave vector and a is the size of 
the molecule. I find how the probability of the molecule 
absorbing an Ej and Mj multipole photon scales with 
ka. In this treatment, I neglect all relativistic interac- 
tions containing electron and nuclear spins. 



II. MOLECULAR AND FIELD HAMILTONIAN 



coordinate systems. The first one transforms as coordi- 
nate components and the second is a geometric quantity 
that does not transform. 

The Hamiltonian for the diatomic molecule can be 
written as [g, |8| 



molecule 



T + V e - e + V n - n + V e -t, 



(2) 



The kinetic energy operator for electrons and nuclei is 
given by 



n p>2 p2 

^ 2m *-< o 2M a 

i—l a— 1,2 



(3) 



where m is the mass of an electron and Mi and M2 are 
the masses of the nuclei, and n is the number of elec- 
trons. Note that I have not assumed that the molecule is 
electrically neutral. The electron and nuclear momentum 
operators, Pi = gp- and P Q = ^p— , are conjugate to the 
coordinates, R; = (Xi,Yi, Zi) and R Q = (X a , Y a , Z a ). 
The Coulomb interaction terms for electrons, nuclei, and 
electrons and nuclei, respectively, are given by 



S IRi-R; 



V e -e — 47reo 2^ jl 
Kj 

T/ _ Zj Z 2 e 1 

v n -n — 47TEQ iRa-Rj 

2 n Z 

o=l,2i=l 



(4) 



where e is the electronic charge, £0 is the permitivity of 
vacuum, Z a is the atomic number and M a is the mass of 
nucleus a. 



I consider a closed system consisting of a diatomic 
molecule interacting with a quantized radiation field. 
The Hamiltonian can be described as a sum of three 
terms [I [Ml: 

H Hmolecule + H field + H int (t) (1) 

where H mo i eC uie is the Hamiltonian for the molecule alone 
(with no field present), H field is the Hamiltonian for the 
field (with no molecule present) and Hi nt is the inter- 
action Hamiltonian between the radiation field and the 
molecule. I assume the diatomic molecule consists of n 
electrons and two nuclei. I use a laboratory system of 
Cartesian coordinates, Ri = (Xj, Yi, Zi), i — 1,2, ... ,n, 
for the electrons and R Q = (X a ,Y a ,Z a ), a — 1,2, 
for the nuclei. I take the basis vectors of the Carte- 
sian coordinate system to be (e x ,e y ,e z ) = (ei, 62,63). 
The position of the origin of the laboratory system 
of coordinates is arbitrary. In what follows, I distin- 
guish between coordinates, e.g., coordinates of the ith 
electron, Rj = (Xi,Yi, Z t ), and the geometric object, 
(Xie x , Yie y , Zie z , ), which represents the position of the 
ith electron and is independent of the coordinate system. 
Note that these two quantities, Ri = (Xi,Yi, Zi) and 
(Xie x , Yie y , Zi& z , ), transform differently under change of 



III. PHOTONS WITH WELL-DEFINED 
ANGULAR MOMENTUM 

For the the radiation field, I can choose Laguerre- 
Gaussian modes [H|,[2(|, which are solutions of the parax- 
ial wave equation, and in which a separation can be made 
between spin angular momentum and orbital angular mo- 
mentum. These types of photons carry angular momen- 
tum hi. Such modes are convenient when dealing with 
cylindrical symmetry, however, they are only approxi- 
mate solutions of Maxwell equations. Alternatively, we 
can choose vector spherical harmonic modes 0, 0] , which 
are exact solutions of Maxwell's equations, but then spin 
angular momentum and orbital angular momentum can- 
not be separated. This is not a flaw of vector spheri- 
cal harmonics, rather it is a general feature of electro- 
magnetic radiation fields. The vector spherical harmonic 
modes are a complete set of vector eigcnfunctions, so any 
field can be expressed in terms of these modes. For these 
reasons, I take the radiation field to be described by vec- 
tor spherical harmonic modes. The Hamiltonian for the 
radiation field, see Eq. (fTJ), is given by 

H field = ^2 hoJk ("ifcjm a\kjm + -) (5) 



3 



where a Xk ^ m and a\kj m are the creation and destruction 
operators for electric or magnetic multipole photons, la- 
beled by A = (E,M), respectively, having wave vector 
magnitude k — uJk/c, and carrying angular momentum 
jh and angular momentum projection mh on z-axis . The 
creation and annihilation operators satisfy the commuta- 
tion relations [3] 



(6) 



Each field state with well-defined photon number has 
well-defined energy. Each photon in a state labeled by 
quantum numbers (A, k,j, m) has wave vector k, square 
of angular momentum j(j + l)h 2 , z-component of an- 
gular momentum mh, and parity (— 1) J for A = E and 
parity (— 1) J+1 for A = M. Photons are of two types: an 
Ej photon is an electric photon and a Mj photon is a 
magnetic photon, each carrying jh units of angular mo- 
mentum [j, H, El ■ The vector potential operator is given 
by 



A(r,t) 



E 

\,k,j,m 



O-Xkjm -&kjm 



(r,t) 



Xkjrn ' 



( 7 ) 

where A^- m (r, t) is the photon wave function for elec- 
tric 2 J pole radiation and Ajy TO (r, t) is the photon wave 
function for magnetic 2 J pole radiation [H, H, H|. In a 
large sphere of radius 1Z, the magnetic and electric pho- 
ton wave functions are given by [4j 



direct products of two different irreducible representa- 
tions of the three-dimensional rotation group, e^, and 
Yim{Q,<fr), coupled by Clebsch- Gordon coeficients: 



Yj^ m (6, ' 



(1,1, n, m - (i | jm) e^F;, 



771 — fl 



(12) 



Here are three spin 5 = 1 eigenfunctions, X^(cr), 
whose components are e M = {x M (+l), X«(0), Xm(-1)}, 
and a = +1,0,-1 is the spin variable @, Q. The e^, 
= {— 1,0, +1}, are the spherical basis vectors, which 
are simultaneous eigenfunctions of the spin angular mo- 
mentum operators S 2 and S z : 



S 2 e, 



2e p 

I- 1 e /i 



(13) 



The spherical basis vectors are related to the Cartesian 
basis vectors, {e x ,e y ,e z }, by 

1 , 

e + = --^=(e x + ie y ) 
e = —j= (e x — ie y ) 



eo 



(14) 



In terms of the Cartesian basis vectors, {e x ,e y ,e z }, the 
action of the spin operator is given by 



Sie k = i ^2 eua e/ 

2=1,2,3 



(15) 



Afcj m (r, t) 



fi(J<T)Y jjm (0,4,)e- 



(8) 



and 



Afcj m (r, t) 



\fj fj+l( kr ) Y 3,j+hm(6,' 



■VT+T/i-rCfcOY^-r, 



k (9) 



where photon frequency is u>k — ck and the spherical 
quantization volume is V. The radial functions, fj(k 
are given in terms of Bessel functions of the first kind 



fi(kr) 



(2V) 



1/2 



1 



J l+ i (k; 



R 



J, 



(10) 



and the allowed set of discrete wave vectors, ki n , where I 
and n are positive integers, are given by 



(h n K)=0 



(11) 



These discrete wave vectors determine the photon ener- 
gies, huJk = hcki n , in a sphere of radius R. For a given 
value of I, the spectrum of allowed wave vectors, {kin} 
is the same for electric and magnetic photons. 

The three vector wave functions (vector spherical har- 
monics), Yjj^ m (9, 0), j = I + 1, 1, \l — 1|, are defined as 



where aki is the totally antisymmetric Levi-Civita ten- 
sor, where ei23 = +1. 

The Yi >m in Eq. (|12l) arc the usual scalar spherical Har- 
monics that are eigenfunctions of the orbital angular mo- 
mentum operators L 2 and L z : 

L 2 Y lm = l(l + l)Y lm 

L z Y lm = mY lm (16) 

According to the definitions in EQ. (fT2|) - (fT6]) . the product 
functions, Yi tJn , are simultaneous eigenfunctions of L 2 , 
L z , S 2 and S z , with eigenvalues 1(1 + 1), m, 2, and \i, 
respectively. 

The spherical vector harmonics, Yjj tm (9,(f>), form an 
irreducible representaion of the three-dimensional rota- 
tion group on the unit sphere, and are simultaneous 
eigenfunctions of J 2 , J z , L 2 , and S 2 : 

(17) 

) (18) 
(19) 
(20) 
4>), satisfy 



JzY j,l,m (9, 


4>) = mYj^m (9, 4>) 


J 2 Yjj }Tn (6, 


4>) = i(i + i)Yi,j, m 


L 2 Yjj^ m (8, 


0) = 1(1 + l)Y jXm 


S Yj t i yTn (6, 


4>) = V£ jthm {6,<l>) 



The vector spherical harmonics, Yjj jT 
the orthogonality relation 



d 2 nY 



jlm 



j'l'rn' 



— Sjj< 5lV $n 



(21) 
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where the integration is over the unit sphere, and com- 
pleteness relation 



- 3,1,171 1 



- j : l,m 1 



= I 5 2 (fi - J!') (22) 



jlm 



Under the parity operator (spatial inversion), the basis 
vectors change sign, and we take the scalar spherical 
harmonics, Yi m , to be defined so that they are multiplied 
by (—1)', and therefore the vector spherical harmonics 
Yj i m (9, 4>) are eigenfunctions of parity with eigenvalues 

From Eq. (jHJ we see that magnetic photons have 
well-defined spin and orbital angular momentum, since 
Ajy m (r, t) is an eigenfunction of spin and orbital angular 
momentum. However, from Eq. (jH]) we see that electric 
photons with quantum numbers (J, m) are a superposi- 
tion of orbital angular momentum I = j ' — 1 and I = j + 1, 
and do not have well-defined orbital angular momentum. 

The number operator for photons of type A is given by 



(23) 



IV. MOLECULE-FIELD INTERACTION 

The electromagnetic field couples to the charges in the 
molecule, which consist of electrons and nuclei. The in- 
teraction Hamiltonian describing the couplin g of the ra- 
diation field to the molecule is given by [g, I9l4l2l [lij 

n 2 

H int (t) = — V A(Rj, t) ■ Pi - e V -fA(R Q ,t) • P Q 

712 i=l a=l Mq 

(24) 

where the vector potential operator is in laboratory co- 
ordinates, and is given in Eq. (J7J. In Eq. (|2~4|) . the first 
term is the coupling of the field to electrons and the sec- 
ond term is the coupling of the field to the nuclei. The 
radiation field, A(r, t), is quantized in the laboratory co- 
ordinates, R. The electron and nuclear momentum oper- 
ators in Eq. (j2~4")l . and P Q , are also in laboratory coor- 
dinates. I have dropped terms that are quadratic in the 
vector potential because they are responsible for photon 
scattering and are small for photon absorption and emis- 
sion. Also, I have dropped small relativistic terms that 
couple the electron and nuclear spins to the field [!, 0, Q ■ 
Furthermore, in Eq. (|24|) . I have assumed the Coulomb 
gauge for the vector potential, so that div A(r, t) = 
and the scalar potential is zero. 

In many experiments, the wavelength of the electro- 
magnetic radiation field is long compared to the dimen- 
sions of the physical system, such as the molecule. We 
can make a Taylor series expansion of the vector potential 
terms A(Rj, t) and A(R a , t) about the center of mass of 
the molecule R G . Keeping only the first term, I make the 
approximations 



A(Ri,t)« A(R Q ,i) w A(Ro,t) 



(25) 



V. TRANSFORMATION TO CENTER OF MASS 
OF THE MOLECULE 

In order to explore the transfer of angular momen- 
tum from photons to molecules, we must look at de- 
grees of freedom that can "absorb" angular momentum. 
To this end, I transform the molecular Hamiltonian in 
Eq. ([2]) from the laboratory coordinates to center-of- 
mass-of-molecule coordinates. I take the Cartesian basis 
vectors in laboratory coordinates to be collinear with the 
basis vectors in center-of-mass-of-molecule coordinates. 
I follow the notation used by Brown and Carrington Q- 
The position of the center of mass of the molecule, in 
laboratory coordinates, is given by 



n ^2 

R ° = m ^ Ri + M ^ MaKa 

i—1 a— 1 



(26) 



where the total molecular mass is given by 



M 



^2 m i + X! Ma = nm + M l + M 2 ( 27 ) 

a=l 



where n is the number of electrons in the molecule. The 
transformation of electron and nuclear coordinates, from 
laboratory coordinates, Ri and R Q , to the center of mass 
of molecule coordinates, R'j and R' Q , is given by 



R i P-o 
PL rv — P*tv R*n 



(28) 



Note that the transformations in Eq. (|28p are transforma- 
tions of vector components, see comment above Eq. ([2"]l, 
I introduce the inter-nuclear vector, R, in terms of the 
coordinates of the two nuclei (in laboratory coordinates) 



R — R9 — Ri 



(29) 



Using the above definitions of the transformation to the 
center-of-mass-of-molecule coordinates, the transforma- 
tion of the kinetic energy of the molecule, given in Eq. ([3]) , 
to the center-of-mass-of-molecule coordinates, gives [7[ 



T : 



M 2u R 2m ^ 4 2M ^ 

i=l j,k=l 



2M 



^ P' P' ( 



where the momentum operators 

h d 



R 



i dH 

h d 



P' =- 9 



(30) 
(31) 

(32) 

(33) 
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are conjugate to the coordinates, R , R and R^, and 
where fi is the reduced nuclear mass 



1 
/' 



1 

Mi 



1 

M2 



(34) 



The first term in Eq. (|3"0"f is the center of mass kinetic en- 
ergy of the molecule, expressed in laboratory coordinates. 
The second term in Eq. (|30[) is the relative kinetic energy 
of the nuclei, expressed in terms of the inter-nuclear vec- 
tor R. The third term in Eq. (|30l) is the kinetic energy of 
the electrons. The last term in Eq. (|3U)) is the so-called 
mass polarization terms and results from fluctuations of 
electron position about the massive nuclei. 

The transformation of the Coulomb interaction terms, 
given in Eq. ((4]), from the laboratory coordinates to 
center-of-mass-of-molecule coordinates, is done similarly. 



VI. TRANSFORMATION TO CENTER OF 
MASS OF THE NUCLEI 

In order to discuss rotation of the molecule about the 
line joining the nuclei, I define the position of the cen- 
ter of mass of the nuclei (in center-of-mass-of-molecule 
coordinates) : 



R 



1 



N 



Mi + M 2 



(AfiR'i + M 2 R' 2 ) (35) 



The positions of the electrons in the center-of-mass-of- 
nuceli coordinates are given by 

R i = Rj — R n 

1 ™ 

= R, '+ - - , «f Z>j r 'j ( 36 ) 



Mi + M 2 



j'=i 



In the center-of-mass-of-nuceli coordinates, the mo- 
mentum operator conjugate to the electron position, R"j, 
is given by 



h d 



i dK"j 



(37) 



The transformation of nuclear coordinates, R Q , to 
center-of-mass-of-nuceli coordinates is similar. 

The total kinetic energy of the molecule, given in 
equation Eq. (|30[) . transformed to center-of- nuclei coor- 
dinates, is given by @, H3] 

i i i ™ 

T = — V 2 + —P 2 R + —Y P"- 

2M ° 2u R 2m ^ 1 

r i=i 



l n 

t p " p' 

2(M 1 +M 2 ).^ i 1 



(38) 



Transforming the Coulomb terms in Eq. @ and (0| 
to the center-of-mass-of-the-nuclei coordinates (which are 



designated by two primed superscripts), these terms be- 
come 



£ r4- 



1/ Z\ Z 2 e 1 

v n -n — 47r£o R 



(39) 



V e -n — 4 ^ J2 mj\ 



where |R"ij| = |R"i — R"j| is the distance between elec- 
trons, |R"j a | = — R"q| is the distance between the 
«th electron and ath nucleus, and R is the inter-nuclear 
vector defined in Eq. (|29]). with magnitude R — R. Note 
that the R"i Q are expressible in terms of the electron co- 
ordinates, R"i, and the inter- nuclear vector, R, so the 
molecular wavefunction depends on the center of mass 
coordinate, R D , the inter- nuclear coordinate, R, and the 
electron coordinates R"j, for i = 1, 2, • • • , n. 

In summary, when the molecule Hamiltonian in labo- 
ratory coordinates, given in Eq. ([2]), is transformed into 
the center-of-mass-of-nuceli coordinates, it is given by the 
sum of Eq. (|3"8"|) and d3U). Note that in the absence of 
the electromagnetic field, the center-of-mass motion of 
the molecule decouples from the internal motion of the 
molecule, i.e., the molecular Hamiltonian does not con- 
tain mixed terms with center-of-mass and internal de- 
grees of freedom. So, in the absence of the radiation 
field, the total molecular wavefunction is the product: 

tpQLM (Ro) <t> (R, R"i, R" 2 , • • • , R" n> ) (40) 
where the center-of-mass wavefunction is of the form 

i>QLM (Ro) = R Q L {QRo) Y LM (0o, to) (41) 

and 4> (R, R"i, R" 2 , • • • , R" n ,) is the wave function for 
internal degrees of freedom. 



VII. INTERACTION HAMILTONIAN IN CM 
OF NUCLEI COORDINATES 



As previously stated, in the absence of the radiation 
field, the molecule is described by a Hamiltonian that is 
the sum of Eq. (|55|) and (f3"9"| . In this case, the center of 
mass motion, described by coordinate R„, decouples from 
the internal motion, which is described by coordinate R. 
However, when an external radiation field is present, see 
the interaction Hamiltonian given in Eq. (|4"12|) . the ra- 
diation field introduces a coupling between the center of 
mass motion of the molecule and its internal motions 

Next, I assume that the molecule interacts with ra- 
diation, according to the interaction Hamiltonian in 
Eq. It is possible to re- write the molecule-field inter- 
action in a form that allows more insight than that given 
by the interaction Hamiltonian in Eq. (l24l) . In order to 
do this, I make the long wavelength approximation, given 
in Eq. (pZS"]). in the interaction Hamiltonian in Eq. (|24l) . 
Then, I transform the electron and nuclear momentum 
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operators to center-of-mass-of-nuclei coordinates. The 
molecule- field interaction Hamiltonian H int (t) in Eq. (|2"4")l 
then becomes 

H- mt {t) = -^{n -Z x - Z 2 ) A(R OJ t) • P + 



i=l 



m Mi + M 2 



A(R G ,i) • P_r 



(42) 



where P is the relative momentum operator for the cen- 
ter of mass of the molecule (in laboratory coordinates) 
that is defined in Eq. (|3"TT) . P"i is the momentum oper- 
ator of the zth electron (in center-of-mass-of-nuclei co- 
ordinates) given in Eq. (1371) . and P^ is the momentum 
operator that is conjugate to the inter-nuclear vector R 
and is given by Eq. (|2l)|) . 

The total Hamiltonian for the molecule, the field, and 
their interaction is given by the sum of Eq. (|3"5]) . (j3"5)) and 
(|42]) . The coordinates in the total Hamiltonian are the 
center-of-mass position (in laboratory coordinates), R G , 
the inter-nuclear position (nuclear relative coordinate) R, 
and the electron positions in the center-of-mass-of-nuclei 
coordinates, R"j. The basis states for the combined sys- 
tem of molecule and field are of the form 

^ (R , R, R"i, . . . , R"„) |M Q1 > 2 ) Q2 ,---) (43) 

where "J (R D , R, R"i, . . . , R" n ) are the molecular states 
that are eigenstates of the sum of Eq. (|3"5|) and (I3T)1) . 
and |(^i) Ql , (^2)^ , • • •) are the photon states that are 
eigenstates of Eq. ([5]), where there are rii photons in mode 
on and I use the short-hand notation a — {A, k, j, m}. 

The first term in the interaction Hamiltonian in 
Eq. (|42p is an ionic term, and is proportional to the pre- 
factor n — Zx — Z<i. This pre- factor is zero for neutral 
molecules, whose electron number is equal to the sum of 
the atomic numbers of the two nuclei, n — Z\ + Z 2 . For 
the case of a molecule that is not neutral (an ion) this 
first term couples the momentum of the center of mass 
of the molecule to the radiation field. This term leads 
to center of mass motion of the (charged) molecule as a 
whole. 



The second term in the interaction Hamiltonian in 
Eq. (l4"2"j) couples the sum of electron momenta (in center- 
of-mass-of-nuclei coordinates) to the radiation field. This 
term allows a transfer of angular momentum from the 
radiation field to the electrons. The pre-factor for this 
second term is never zero. 

The third term in the interaction Hamiltonian in 
Eq. ((32) couples the momentum, 



Pr = -in 



OR 



x 



8Ry 



dR, 



(44) 



which is conjugate to the inter-nuclear vector, 
R = (Rx, Ry, Rz), to the electromagnetic radiation 
field, A(R D , t). Here, (e x ,e y , e z ) are the basis vectors for 
the space- fixed coordinates of R = (R x ,Ry,Rz)- The 
pre-factor for this term is zero for homo-nuclear diatomic 
molecules which have Mi = M 2 and Z\ = Z 2 . 

Homo-nuclear diatomic molecules have an inversion 
symmetry about the nuclear center of mass. Molecules 
which have Z\ = Z 2 but Mi 7^ M 2 , (e.g., different 
isotopes, such as HD), lack this inversion symmetry, and 
then this third term in Eq. (|42[) is generally non-zero. 
This third term contains a coupling of the radiation field 
to the rotation of the molecule about the line joining 
the nuclei. This can be seen by considering the unitary 
transformation, Ai, from space- fixed coordinates with 
origin at the center-of- nuclear mass, X = (X, Y, Z), 
to molecule-fixed coordinates, x = (x, y, z), where the 
z-axis is pointing along the inter-molecular vector 
R 2 Ri: 



X 

Y \ =M{<t>,6, X ) 
Z 



(45) 



where (0, 9, x) are the Euler angles between coordinates 
X = (X, Y, Z) and x = (x, y, z). I take the sequence of 
rotations through angles {</>, 0,x} to be along the z, y, z 
axes, so the rotation matrix is given by @, H3, [H| 



cos 4> cos 9 cos x ~ sin 4> sin x ~~ sin 4> cos x ~ cos 4> cos 9 sin x cos cf> sin ( 
M.(cj>, 9, x) = sin 4> cos 9 cos x + cos 4> sin x cos cf> cos x — sin <j> cos 6 sin x sin 4> sin t 
— sin 9 cos x sin 9 sin x cos 9 



(46) 



For internuclear vector R = (Rx,Ry,Rz), the trans- coordinates, (Rx, Ry,Rz) — > (R, 9, <j>, x), is given by 
formation from space-fixed coordinates to molecule-fixed 



Rx 

Ry I =M(#AX) 
Rz 





R 



(47) 
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For a diatomic molecule, since the vector 
R = (Rx, Ry, Rz) is parallel to the molecule-fixed 
z-axis, the components (Rx,Ry-,Rz) do not depend 
on x, and we are free to choose its value for conve- 
nience @, USUI- Using the transformation in Eq. (|47|. 
the relative momemtum in Eq. (|44|) becomes |7|: 



Pr = S 



1 ™ ■ d 
n x N — m— — 

R OR 



(48) 



where n is the inter-nuclear unit vector, given in space- 
fixed coordinates by 



R2 R-i 
|R 2 -Ril 



sin 9 cos </> e x + sin 9 sin <f> e y + cos 9 e z 

(49) 

and N is the angular momentum operator in space-fixed 
coordinates with origin at center-of-mass of the nuclei 



N 



cot 9 cos i 



cot 9 sin . 



b d(f> 
,— 



cos . 



— 



c sec 9 cos 



c sec sin < 



d ' 

ox. 



d 



(50) 



Using the form of Pr in Eq. 
Eq. (@2J, gives 

A(R )- P R ,t = -A(R ,t)- 



in the third term in 



1 • 9 

-nxN + J „— 



(51) 



which shows that the radiation field couples to both 
molecular rotations and molecular vibrations. The first 
term in Eq. (|51[) is a coupling of the radiation field to 
molecular rotations, through the molecular angular mo- 
mentum N. The second term in Eq. (|5ip is a coupling 
of the radiation field to molecular vibrations through the 
length of the inter-nuclear distance R. 

Babiker et al. investigated the interaction of a hydro- 
genic bound system consisting of a nucleus and a bound 
electron, interacting with a linearly polarized radiation 
field carrying orbital angular momentum (29| . They 
found that an exchange of orbital angular momentum 
in an electric dipole transition occurs only between light 
and the center of mass motion of the atom — they con- 
cluded that the internal electronic-type motion does not 
participate in any exchange of orbital angular momen- 
tum in a dipole transition. The presence of the second 
term in Eq. (l4"2l suggests that this is not the case for a 
diatomic molecule. 

Alexandrescu et al. studied a charged diatomic 
molecule (an ion, with two nuclei) interacting with a 
Laguerre-Gaussian beam carrying orbital angular mo- 
mentum [30| | . They seem to disagree with the conclusions 
of Babiker et al., stating that the electronic motion and 
the electromagnetic field can exchange one unit of OAM 
within the electronic dipole interaction. 

I note that the interaction Hamiltonian in Eq. (|4"2")l is 
not of the form d • E, which is a dipole coupling to 



the radiation field E. The interaction Hamiltonian in 
Eq. (l4"2"j) contains all multipole couplings with the radia- 
tion field. The interaction Hamiltonian in Eq. (|4"2l allows 
investigation of transfer of angular momentum from the 
radiation field to the molecule, including transfer of angu- 
lar momentum to the rotation of the nuclei. Regarding 
the dipole moment, consider the dipole moment of the 
molecule, expressed in terms of variables R"j and R" Q 
in the center-of-mass-of-nuclei coordinates: 



(52) 



a=l 



When the dipole moment is transformed to the set of 
coordinates {R G , R, R"i, ■ ■ ■ R"n}, it becomes 



1 + I7 (Zl + Z2) l^ R 



+ e(Z 1 + Z 2 )R - e (i|-:|) M R (53) 

In order to isolate the rotational part of the Hamilto- 
nian, contained in Pr, it was necessary to transform to 
coordinates {R Q , R, R"i, . . . R"„}. However, after per- 
forming his transformation, the Hamiltonian in Eq. (|42[) 
is not of the form of a dipole coupling to the radiation 
field, — d • E. This is not a result of the transformation 
itself, rather it is due to the added rotational degrees of 
freedom of a molecule as compared to an atom. 



VIII. TRANSITION PROBABILITIES 

The form of the molecule-radiation interaction in 
Eq. ((32) shows that there is a clear coupling of the ra- 
diation field to molecular rotations. The center of mass 
coordinate, R Q , plays a central role in this coupling. As 
stated earlier, in order to determine the effect of this 
coupling on observable quantities, transition matrix ele- 
ments and probabilities for photon absorption must be 
calculated. In order to estimate these matrix elements, 
we must have knowledge of the form of the wavefunctions 
for the diatomic molecule, which are schematically shown 
in Eq. (|43|) . We must also assume some form for the 
modes of the radiation field, for example, vector spherical 
harmonics or Laguerre-Gaussian modes. Group theoretic 
methods based on tensor operators can then be used to 
determine ratios of transition matrix elements @, H3, HH ■ 
This work is in progress. However, some statements can 
be made about parity selection rules and how the transi- 
tion matrix elements scale. 

Exact selection rules based on conservation of parity 
can be written down. Assume that under inversion of 
coordinates (under the parity operator), the parity of the 
initial and final molecular wavefunctions are Pi and Pf, 
respectively, and the parity of the emitted photon is P 1 , 
where parity eigenvalues takes values ±1. If the initial 
state of the whole system has one photon, and the final 
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state has no photons, then parity conservation requires 
Pi P 1 = Pf, or, since each parity eigenvalue is ±1, parity 
conservation requires Q 

Pi Pf = Pj (54) 

Electric Ej photons have parity P®i = (— 1) J and mag- 
netic Mj photons have parity P* /j = (— For ab- 
sorbing one photon of a given type, parity conservation 
requires 



PiPf = (~l) 



electric Ej photon 



(55) 



or 



PPf = (-l) 



3+1 



magnetic Mj photon (56) 



An estimate of the scaling of transition probabilities 
can be made by considering the molecule as a system 
of charges interacting with a radiation field. Consider a 
molecule in an initial state with energy Ei, absorbing a 
single photon of energy huk = hck, and making a tran- 
sition to a final state with energy Ef. The transition 
matrix element contains the quantity kr, which enters 
in the vector potential in Eqs. ([8} and (|9]). The quan- 
tity r has a characteristic scale of a diatomic molecule, 



on the order of a 



10 



-10 



m. For transitions be- 



tween rotational states, I take the frequency / ~ 300 
GHz. Therefore, there is a small dimensionless quantity, 
kr ~ ka 2ira/X ~ 6 x 10~ 7 in the transition matrix 
elements. The vector potentials in Eqs. © and © de- 
pend on the function fi(kr), which in turn depend on 
the Bessel function J l+ i(kr). Using the small argument, 
kr -C 1, expansion of the Bessel function, in Eq. (|TU| . 
the radial function, fi(kr), scales as (up to a constant): 



fi(kr)~(kr) 



1 



1 • 3 • 5 - - • (2Z + 1) 



(57) 



The vector potential for Ej photons scales as 
A^ m (r) ~ fj_x (kr), and the vector potential for Mj 
photons scales as Aj!j TO (r) ~ fj (kr). The probability 
of absorbing a photon scales as the square of the vec- 
tor potential (square of the matrix element). Therefore, 
the ratio of the probability for absorbing an Mj photon, 
P(Mj), divided by the probability of absorbing an Ej 
photon, P(Ej), scales as 



P(Mj) (kaY 



4 x 10 



-13 



P(Ej) (2j + l) 2 (2j + l) 



(58) 



For a given angular momentum j, the probability of ab- 
sorption of a magnetic photon is much less than absorp- 
tion of an electric photon. 

The ratio of the probability for absorbing an E(j + 1) 
photon, P(E(j + 1), divided by the probability of absorb- 
ing an Ej photon, P(Ej), also scales as 



P (E(j + 1)) (kaf 4xl0" 13 
P (Ej) (2j + 1) 2 ~ (2j + l) 2 



(59) 



Form Eq. (I59|) . we see that the probability of absorb- 
ing higher angular momentum photons rapidly decreases 
with increasing j. As the simplest example, the ratio of 
probability of absorbing an E2 photon to the probability 
of absorbing an El photon is: 



P(E2) 
P(E1) 



4 x 10" 



(60) 



The scaling of the absorption probabilities for 
molecules, given in Eqs. (|58p - (j6"0")l . makes it a challeng- 
ing proposition to observe absorption of photons having 
angular momentum j > 1. 



IX. SUMMARY 

I have investigated the interaction between a diatomic 
molecule and photons with well-defined angular momen- 
tum. I have exploited the transformation from labora- 
tory coordinates to coordinates with origin at the center 
of mass of the nuclei to show the nuclear rotational de- 
grees of freedom explicitly (through the vector R). For 
the photon modes, I have assumed spherical vector har- 
monics in a large sphere of radius R. These modes are 
exact solutions of Maxwell equations, and consequently, 
any other modes, such as the Laguerre-Gaussian modes, 
can be expanded in terms of them. However, in general, 
spin and orbital angular momentum are not separable in 
vector spherical harmonic modes. I used the standard 
minimal coupling between elementary charges and radi- 
ation field to couple the radiation field to the molecule, 
neglecting quadratic terms in the vector potential and ne- 
glecting relativistic interaction with spins. For molecules, 
which absorb at millimeter wavelengths, there is a small 
parameter, k a where k is a wave vector and a is the size 
of the molecule. This small parameter was used to deter- 
mine the scaling of absorption probabilities of Ej and Mj 
type photons. For transitions between rotational levels of 
the molecule having AJ > 1, the probability of absorp- 
tion of photons with angular momentum j > 1 is small, 
because the factor fca< 1, where wave vector k = 2ir/\ 
corresponds to photons in the millimeter wavelength re- 
gion, see Section IVIIII Note that the small parameter, 
ka, also exists for atoms [l9l - l2l| . however, it is larger 
for atoms because k corresponds to visible wavelengths, 
whereas for rotational transitions in molecules k corre- 
sponds to millimeter wavelengths. 

For electric photons having arbitrary angular momen- 
tum jh, the probability of rotational transition by ab- 
sorbing an electric E(j + 1) photon divided by the prob- 
ability of absorbing an electric Ej photon, scales as 
(ka) 2 /(2j + l) 2 . The probability of absorbing a mag- 
netic Mj photon is much smaller than the probability 
of absorbing an electric Ej photon. The probability of 
absorbing a magnetic Mj photon, divided by the prob- 
ability of absorbing an electric Ej photon, also scales as 

M7(2j + i) 2 . 
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